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For a long time, it is believed that the light by light scattering is described properly by the 
Lagrangian density obtained by Heisenberg and Euler. Here, we present a new calculation which 
is based on the modern field theory technique. It is found that the light-light scattering is com- 
pletely different from the old expression. The reason is basically due to the unphysical condi- 
tion (gauge condition) which was employed by the QED calcualtion of Karplus and Neumann. 
The correct cross section of light-light scattering at low energy of (— <C 1) can be written as 

S^(6^d^P + ^«'"''' + -''"- 
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I. INTRODUCTION 

It has been believed that photon-photon scattering can be described by the Lagrangian density of Heisenbergand 
Euler [H • In addition, this cross section is confirmed by the QED calculation of Karplus and Neumann @1 ■ 
However, the calculation by Karplus and Neumann is not reliable since they have put some additional conditions on 
the QED calculation so as to reproduce the result obtained by Heisenberg and Euler. 

In this respect, it is important that the proper treatment of the photon-photon scattering should be made again at 
the present stage. The calculation of photon-photon scattering itself is straightforward since the Feynman diagrams 
of the fourth order perturbation calculation can be done without any difficulties. The interaction Lagrangian density 
of H' = —ejfj,A'^ is gauge invariant as long as the fermion current is conserved, which is indeed the case. Therefore, 
the QED calculation itself has no conceptual difficulty and therefore we should carry out the S-matrix evaluation of 
photon-photon scattering. 
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II. FEYNMAN AMPLITUDE OF PHOTON-PHOTON SCATTERING 

The Feynman amplitude of the photon-photon scattering can be written as 

^rr'ss' ^ 2{Ma + Mb + Mc) (1.1) 

where we note that the Feynman amphtude is not affected by the direction of the loop momentum Q • The amplitude 
Ma can be explicitly written as 



(27r)4 J ^ [{q - I - k)'^ - m'^ + ie] [{q - I'Y - m? + is] [q^ - + is] [{q - l)^ - m'^ + ie] 

xel{k)4{k')ei{l')el{l). (1.2) 

Here, it should be important to note that the total Feynman amplitude of the photon-photon scattering does not 
have any logarithmic divergence even though the amplitude Ma alone has the logarithmic divergence. This is quite 
important in that the physical processes should not have any divergences, and indeed the photon-photon scattering is 
just the case. This strongly suggests that the evaluation of the Feynman amplitude of the photon-photon scattering 
process must be directly connected to the real physical process which should be observed by experiments. 

HI. DIVERGENCE 

The amplitude Ma has the logarithmic divergence, and this can be seen when we check the large q behavior [4]. In 
this case, we find from (1.2) 

^''^~J2^J "^^J^^^"^^'" + ' ^a-'d'^Kelei el (2.1) 

which has obviously the logarithmic divergnce. However, if we add all of the amplitude together, then we find from 
(1.1) 

+5^V^+5^^^?"^"-2<7■^V" (2.2) 

= 0. 

This means that the total amplitude has no divergence at all because of the cancellation, and it is indeed finite. 
Therefore, we do not have to employ any specific regularization scheme, and thus the evaluation is very reliable. 

IV. DEFINITION OF POLARIZATION VECTOR 

Here, we take the polarization vector as defined by Lifshitz Q 

f(l) - f(l) - f(l) - f(l) - ^ ^ .(2) - ^ a X F^l)^ - 

^4 ^ — ^4 "^j — ^3 
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denotes the polarization vector of photons. Here, we take the Coulomb gauge fixing V • ^1 = 0. 
Each photon has the following momenta 

initial state 

photon 1 : fc^ = (w, k) photon 2 : = (w, I) 

final state 

photon 3 : I'f" ^ {uj,l') photon 4 : fc'^ = (w,fc')- 

uj=\k\^ \l\ = Ife'l = 1^1 
Also, by noting that there is no rest system, we find 

l^-k I' = k'. 

V. CALCULATION OF Ma AT LOW ENERGY 

Now, we carry out the calculation of Ala as an example 

Ma = Mr'"e;(fc)e:(06i'(/')e;'(fc') (4.1) 

where M^""*"^ can be written as 



(27r)4 7 [(q - I - ky - m'^ + ie][{q - l')-^ - m'^ + ie][q^ - + ie][iq - ly - + le]' 
By making use of the Feynman parameter 8] , we find 

where 

Af" = l^'zl + (/' - l)^'z2 + k'^'zs 
B = 2l.kz3 - 



By putting q — A — t, we find 
(27r)4 



= -7H:t3! / d^i / dz2 / dz: 



./O "'0 
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, ^4^ Tr[7^(/ + 4 - f - ^ + m)-f''{lf + + m)-f^{f + 4 + m)Y [j + / - f + m)] 

' [i2 _ + B + ze]4 



Here, for simplicity, we define 
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and thus 



"'0 



21 /•Z2 



3! / dzi / dz2 / dz3 / dH 



[<2 - ^2 + 5 + 



In the numerator, we find Tr [odd- numbers of 7-matrices] — and by noting J dH t^t'^ ■ ■ ■ t\ f{t^) = we find 



odd 



Tr[7^(|^ + i + vi)-i''{if + ^ + m)7^(/ + ^ + 771)7'' + ^ + 777)] 



where 



Now, we find 



__Tr[{7'^7''7^,z!7''^ + -i'' (jL + 7^7''^7^^7'' 

(^Mi^A^ = Tr[7^(^ + 777)7'"(^ + 777)7-^(^ + 777)7''(^+ 777)]. 



(2^)4 



Jo 



3! / dzi / dzi / dz3 / d i 



(^2 - + S + ief 



and we carry out the integration and obtain 

\4 r\ 



(27r) 



Jo Jo 



\B~~A? ^ {B - ^2)2 



g''V'^-2g-V^) In 



A2 



S-yl2 



11 



where A denotes the cutoff momentum. 

Here, we consider the case of — <C 1 and expand the integration in terms of — powers. This approximation is quite 
reasonable as we show it later [J. 

Now, we consider the expansion in terms of ^. Noting that the 771 never appears in fc, I, k' , I' in F^^'^^" , Qi^f^^o- 
but only w appears, we can rewrite F*"^Acr^ Q^ivXa obtain 



where 



uX(7 



^Tr[{7^7'^7^^7''^ + 7^7^^7'^7''^ + 7''7'^^7^^7'' 



+ (7^^7")(7\"^ + 7''^7^ + ^7V)} 



g^^uxa ^ Xr[7^^7"^7V7'^^ 



Here, Q^''-*'"', i?^''-^'^, S"^"-*"^, yMfAa ^^j^g dimensionless. In this case, we can easily expand 

^^uXa j^^yXa g^vXcr rp^yXa 
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Now, we can write 



B 



2:2(2:2 - -21 - zs) sin^ - - Z3Z1 cos^ 2 



and therefore we have 



^0 



Zl /•Z2 

'2 
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Here, we omit the index of /x, z/, A, ct in Q'^"^'^, ii'"'^'^, S'**''^'^, T'*'^^'^. Therefore, we find 



i/Acr 



(2 
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.^2U (l-^2{7)2 
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^In 
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Expanding in terms of ^ = ^ , we find 



>UI>Acr „■ 2 



/ d^i /" ' d22 r dzz[^Q + {R + 8QU)^^ + {RU + SQU^)^'^ + ■ ■ 
Jo Jo Jo I 

+T+{R + 2T[/)^2 + (5 + 2RU + 3TC/2)^4 ■ 



"'0 Jo 



(27r)4 

ient c 

of Ma, Mb, Me- Therefore, we find 
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The coefficient of Q at the lowest order can be written as 8 + 8(ln|^| — which does not depend on the shape 



^r"'^ = -Mli^' f dz, r fj dzs T = T 



(27r) JO ^0 
and we can write it explicitly as 



{2t:Y 6 
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VI. TOTAL AMPLITUDE 

In this way, we can obtain the shape of M^""^" , Mj^"^" and M^"^" at ^ < 1. Thus, we can write for M^, Mb, 



as 

^4 



Mb = + .9'"^.9'^^ - g'-'^g^nyyAe'; (5.1.6) 



Thus, the total amphtude can be calculated as 



g^^g"" + g^^g"^ ~ g'^'g^" 



g'^^g^'^ + g^-g^^ - g-^g^- \eleteie:; 



(2^)4 6 ^ ' ^ ^ ■ ^ ^ ;> 

Therefore, we can write M as 

4 



M = -z - (5^-/'' + gf^-g-^ + g^^g^^) ^d^e^ (5.2) 



where a — 



VII. PHOTON-PHOTON SCATTERING 

Now, the photon-photon scattering cross section can be written as Q 

da 1 1 



MY (6.1) 



dVL 647r2 (2a;)2 

where we sum up the final states of the photon polarization state and make average of the initial polarization states. 



rr' ss' 1 2 

J \1VI 
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rr' ss' 



Among the above, the nonvanishing term is 



1221|2 I |j,f2112|2 I 17.^-121212 



\M''''Y + \M 



|M1212|2 + |m2121|2| 
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Since we know that 



we obtain 



M"22 ^ ^,^2211 ^1221 ^ ^^2112 ^,^1212 ^ ^2121 



i J2 |M"''^^'P = i {|M""|2 + |M2222|2^2|M"22|2^2|Ml221|2+2|Afl212|2l 



rr' ss 



Further, we find 



O 



^j1221 ^ 5;^^2^^^^ 



and thusA 



Therefore, for — ^ 1, we find 



i ^ |M|2 = {3 + 2 cos2 + cos^ 6} . 



'^'^ ^ (3 + 2cos2 + cos4 0) (6.2) 



(67r)2 (2a;) 

2 

where a = . This is the scattering cross section in ^ 1. This result is very different from the one obtained by 
Euler-Heisenberg The calculated result of Euler-Heisenberg can be written as 

da liQa'^ 



(3 + cos2 0)2. (6.3) 

\mJ 



dn (1807r)2m2 Km 

In particular, the energy dependence of the cross section is completely different from each other. At the low energy 
limit, the new cross section becomes larger while the Euler-Heisenberg cross section becomes smaller. At w ~ 1 eV, 
the Euler-Heisenberg cross section with 90 degree becomes 



da 139a^ 



dfl (607r)2m2 \ m 
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9.3 X 10^^^ cm^ ~ 9.3 x 10"'*^ b 



which is extremely small. This suggests that the Euler-Heisenberg cross section is practically impossible to measure 
since it is too small. On the other hand, the present calculation predicts the cross section at w ~ 1 eV 

da 3a^ n i r, 

if; = no 2 - 2-3 X 10-21 cm^ ~ 2.3 x 10^ b. 
dit (127r)^a;^ 

This is rather large, and as far as the magnitude of the cross section is concerned, it should be indeed measurable. 
These problems are discussed in detail in [9|. 
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VIII. WHY IS EULER-HEISENBERG LAGRANGIAN INCORRECT ? 

Here, we briefly discribe the physical reason why the Euler-Heisenberg is incorrect. 



A. Euler-Heisenberg Lagrangian density 



The Euler-Heisenberg Lagrangian density is given as 

2a2 



45m4 



7{E-Bf 



(7.1.a) 



They first calculate the vaccum polarization effects when they apply the electromagnetic fields to the Dirac vaccum 
states which are composed out of negative energy fermions. In this case, they obtained the effects of the vaccum 
polarization, and constructed the effective Lagrangian density which can simulate the results. The basic problem is 
that they treated the vacuum state as if it were dielectric matter. 



B. Karplus and Neuman's calculation 



Now, we discuss the result calculated by Karplus and Neuman who started their calculation from the QED La- 
grangian density. Here, we want to clarify why Karplus and Neuman obtained the same results as that of Euler- 
Heisenberg, in spite of the fact that Karplus and Neuman employed the modern field theory terminology. First, we 
rewrite the Euler-Heisenberg Lagrangian density as 

C = -Y^{5(F,.(x)F^^(a:))^ - UF,,{x)F'^\x)F^,{x)F''^ {x)] (7.1.6) 

where 

F^'^ix) = d^'A^ix) - d'^A^'ix). 

Karplus and Neuman calculated first the QED evaluation based on the QED Lagrangian density. However, they 
always imposed the conditions that they should be able to reproduce the results of the Euler-Heisenberg calculation. 
This must be a mistery why they believed that they should get the same result as the one by Euler-Heisenberg. 
However, this may well be conected to the additional conditions which are often called "gauge condition" even though 
there is no physical reason for this gauge condition (lo| . 



IX. CONCLUSION 



We have presented the new QED calculation of photon-photon scattering at low energy. The calculation is straight- 
forward since there is no logarithmic divergence in the evaluation of the photon-photon scattering Feynman diagrams. 
The result is very different from the Euler-Heisenberg calculation, and therefore the photon-photon cross section 
should be measured by experiments in order to clarify which of the cross section should be preferred by nature. 

X. ACKNOWLEDGEMENTS 

The author is grateful to Prof. T. Fujita for his helpful discussions and comments. 



[1] H. Euler and B. Kockel, Naturwissensch. 23, 246 (1935) 

[2] H. Euler, Ann. der. Physik. 26, 398 (1936) 

[3] W. Heisenberg and H. Euler, Z.Phys. 98, 714 (1936) 

[4] V.B. Berestetskii, E.M. Lifshitz and L.P. Pitaevskii, " Relativistic Quantum Theory", (Pergamon Press, 1974) 



9 



[5] R. Karplus and M. Ncuman, Phys. Rev. 80, 380 (1950) 

[6] R. Karplus and M. Ncuman, Phys. Rev. 83, 776 (1951) 

[7] K. Nishijima, "Fields and Particles", (W.A. Benjamin, INC, 1969) 

[8] F. Mandl and G. Shaw, "Quantum Field Theory", (John Wiley & Sons, 1993) 

[9] T. Fujita and N. Kanda, A Proposal to Measure Photon-Photon Scattering (to be published) 

[10] T. Fujita and N. Kanda, "Tomonaga's Conjecture on Photon Self-Energy", physics.gen-ph/1102.2974 



